package com.mdd.common.utils;

public class RangeCalculationUtil {

	/**
	 * 地球平均半径（单位：米）
	 */
	private static final double EARTH_AVG_RADIUS = 6371000;

	/**
	 * 方法一：（高德地图计算方法）
	 *
	 * @param longitude1 第一点的经度
	 * @param latitude1  第一点的纬度
	 * @param longitude2 第二点的经度
	 * @param latitude2  第二点的纬度
	 * @return 返回的距离，单位m
	 */
	public static Double getDistanceByGaoDe(double longitude1, double latitude1, double longitude2, double latitude2) {
		if (longitude1 == 0 || latitude1 == 0 || latitude2 == 0 || longitude2 == 0) {
			return -1.0;
		}
		longitude1 *= 0.01745329251994329;
		latitude1 *= 0.01745329251994329;
		longitude2 *= 0.01745329251994329;
		latitude2 *= 0.01745329251994329;
		double var1 = Math.sin(longitude1);
		double var2 = Math.sin(latitude1);
		double var3 = Math.cos(longitude1);
		double var4 = Math.cos(latitude1);
		double var5 = Math.sin(longitude2);
		double var6 = Math.sin(latitude2);
		double var7 = Math.cos(longitude2);
		double var8 = Math.cos(latitude2);
		double[] var10 = new double[3];
		double[] var20 = new double[3];
		var10[0] = var4 * var3;
		var10[1] = var4 * var1;
		var10[2] = var2;
		var20[0] = var8 * var7;
		var20[1] = var8 * var5;
		var20[2] = var6;

		return Math.asin(Math.sqrt((var10[0] - var20[0]) * (var10[0] - var20[0]) + (var10[1] - var20[1]) * (var10[1] - var20[1]) + (var10[2] - var20[2]) * (var10[2] - var20[2])) / 2.0) * 1.27420015798544E7;
		// 结果四舍五入 保留2位小数
		//return new BigDecimal(distance).setScale(2, RoundingMode.HALF_UP).doubleValue();
	}


	private static double rad(double d) {
		return d * Math.PI / 180.0;
	}

	/**
	 * 方法二：（基于googleMap中的算法得到两经纬度之间的距离,计算精度与谷歌地图的距离精度差不多。）
	 *
	 * @param longitude1 第一点的经度
	 * @param latitude1  第一点的纬度
	 * @param longitude2 第二点的经度
	 * @param latitude2  第二点的纬度
	 * @return 返回的距离，单位m
	 */
	public static double getDistanceByGoogle(double longitude1, double latitude1, double longitude2, double latitude2) {
		double radLat1 = rad(latitude1);
		double radLat2 = rad(latitude2);
		double a = radLat1 - radLat2;
		double b = rad(longitude1) - rad(longitude2);
		double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a / 2), 2) + Math.cos(radLat1) * Math.cos(radLat2) * Math.pow(Math.sin(b / 2), 2)));
		s = s * EARTH_AVG_RADIUS;
		s = Math.round(s * 10000d) / 10000d;
		return s;
	}
}
